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We investigate the dynamics under diabatic expansions/compressions of linear ion chains. 

Combining a dynamical normal-mode harmonic approximation with the invariant-based inverse¬ 
engineering technique, we design protocols that minimize the final motional excitation of the ions. 

This can substantially reduce the transition time between high and low trap-frequency operations, 
potentially contributing to the development of scalable quantum information processing. 
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I. INTRODUCTION 

Trapped-ion chains provide currently a leading archi¬ 
tecture to test and develop quantum information process¬ 
ing protocols [H-Q . A toolbox of fundamental ion manip¬ 
ulations is required to implement the logical operations 
or simulations. This comprises procedures to control the 
ions’ internal (electronic) and external (motional) degrees 
of freedom, and also to couple them together. Such pro¬ 
cedures should be not only feasible but also fast. Achiev¬ 
ing high speeds is important in itself, to allow for more 
operations per unit time, but it is also instrumental in 
suppressing the effects of noise and decoherence. 

Regarding motional control, a number of operations 
have been identified as relevant for different scalable 
trapped-ion quantum information architectures j in¬ 
cluding ion transport |7l-ll2ll . and ion-chain splitting and 
recombination [l3l4l^ . These operations can be per¬ 
formed on single or mixed-species ion chains [l^ . allow¬ 
ing for sympathetic ion cooling or quantum-logic spec¬ 
troscopy [13 • 

The physical operation that we consider here is fast 
control of the motional frequencies of the trapped ions, 
which in the case of multiple ions leads to chain expan¬ 
sions and compressions. Several elementary protocols 
benefit from a high trap-frequency whereas others are 
better performed with low frequencies. Therefore, a fast 
transition between them without inducing final excita¬ 
tions is a worthwhile goal. 

Operations that benefit from high motional frequencies 
(ie. large potential curvature, small inter-ion distance, 
and small Lamb-Dicke parameters) include: 

• Doppler laser cooling, since the mean phonon num¬ 
ber is lower for tighter traps [l^ : 

• any operation where a single motional normal mode 
(NM) of an ion chain needs to be spectrally re¬ 
solved, since the NM frequency splitting is propor¬ 
tional to the trap curvature [l3|; 

• operations which make use of motional sidebands 
and whose fidelity is limited by off-resonantly driv¬ 
ing carrier transitions on the qubit. 


On the other hand, operations where a lower motional 
frequency is desired include: 

• single-ion addressing in a multi-ion crystal; 


• resolved sideband cooling, which cools at a rate pro¬ 
portional to the square of the Lamb-Dicke param¬ 
eter [l^ : 

• geometric phase gates [2^, which are faster for 
larger Lamb-Dicke parameters. 



In many cases a compromise will be optimal, depending 
on the dominant limitations for a particular experiment. 

Fast expansions/compressions without final excitation 
have been designed in a number of different ways 
Invariant-based engineering or scaling methods 
were realized experimentally for a non-interacting cold- 
atom cloud [1^ and a Bose-Einstein condensate [2^ [2^ . 
However, the methods used rely on single-particles, BEG 
dynamics, or equal masses, and are not directly applica¬ 
ble to an arbitrary interacting ion chain. We propose here 
a method to design trap expansions and compressions 
faster than adiabatically and without final motional ex¬ 
citation. Specifically we define dynamical normal modes 
similar to the ones defined for shuttling ion chains in Q 
and apply invariant-based inverse-engineering techniques 
by either exact or approximate methods. 

We first discuss two-ion chains in Sec. im both for ions 
of equal mass, and ions of different mass, and then extend 
the analysis in Sec. IHII to longer chains. In the examples 
only expansions of the trapping potential are considered, 
as compressions may be performed with the same time- 
evolution of the spring constant, only time-reversed. 


II. 2-ION CHAIN EXPANSION 


We will deal with a one-dimensional trap containing 
an A-ion chain whose Hamiltonian in terms of the posi¬ 
tions {qi} and momenta {pi} of the ions in the laboratory 
frame is 
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where Cc = 4 ^^, with eo the vacuum permittivity. 
uo{t) is the common (time-dependent) spring constant 
that defines the oscillation frequencies U!j{t)/{2TT) for 
the different ions in the absence of Coulomb coupling: 
uo{t) = miUjKt) = 7712^2(t) = • • • = m]s[U!%{t). All ions 
are assumed to have the same charge e, and be ordered as 
qi > q 2 > ■ ■ ■ > qN, with negligible overlap of probability 
densities as a result of the Coulomb repulsion. The po¬ 
tential term V{qi,q 2 ) = H — in the Hamilto¬ 

nian o for two ions is minimal at the equilibrium points 
q^i^ = xq/2, = -a:o/ 2 , where xq = a:o(t) = 

is the equilibrium distance between the two ions. Instan¬ 
taneous, mass-weighted, NM coordinates are defined by 
diagonalizing the matrix 

The time-dependent eigenvalues are [s^ 

^ ^ + 

where we have relabeled mi —^ m and m 2 —>■ ^m, and 
omitted the explicit time dependences to avoid a cum¬ 
bersome notation, i.e., A± = A±(t) and wi = wi(t). The 
time-dependent angular frequencies for each mode are 
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and the eigenvectors corresponding to these eigenvalues 
are v± = (a±,b±)'^, where 



The instantaneous, dynamical normal-mode (mass 
weighted) coordinates are finally 



The quantum dynamics of a state |'0) governed by H in 
the laboratory frame may be transformed into the moving 
frame of NM coordinates by the unitary operator 


U = 


dq+dq_d(7id(?2|q+,q-)(q+,q-ki,92)(<?i,<72|, (6) 


where (q+, q_|gi, 92 ) = S[qi - gi(q+, q_)] 5 [( 3'2 - 

( 72 (q-i-, P-)]- The Hamiltonian in the dynamical equation 
for 1 ^') = 17 1 ^) is given by 

H' = UHU"' - ihU{dtU'') = 

== - PouPu + , (7) 

L/ ^ ■' 

where cubic and higher order terms in the coordinates 
have been neglected, = ±, p± are (mass weighted) 
momenta conjugate to q±, and 


Po± = -q± = 3 (““± 

are functions of time with the same dimensions as the 
mass weighted momenta. They appear because of the 
time dependence of the NM coordinates through xq, 
which is a function of uji{t). These po± functions act 
as momentum shifts in a further unitary transformation 
which suppresses the terms linear in p±, 

— g-i(po+q++po-q-)/fi- 

m = uw), 

H" = UH'U^ - imidtU"^) = 



This Hamiltonian corresponds to two effective harmonic 
oscillators with time-dependent frequencies and a time- 
dependent moving center. Note that the “motion” of the 
harmonic oscillators is in the normal-mode-coordinate 
space, and that the actual center of the external trap 
in the laboratory frame is fixed. According to Eqs. 
m and ([5]) both the NM harmonic oscillators’ centers 
(~Po±/^ 2 ±) and the frequencies (H±) depend on 
This is important as, to solve the dynamics for given 
the oscillators are effectively independent. How¬ 
ever, from an inverse-engineering perspective, their time- 
dependent parameters cannot be designed independently. 
This “coupling” is here more involved than for the trans¬ 
port of two ions in a rigidly moving harmonic trap Q, 
where Po±(I) take different forms which depend on the 
trap position but not on the trap frequency. A different 
approach is thus required. 

The Lewis-Riesenfeld invariants [3l| of the two oscilla¬ 
tors are 


/mi -I- b± y/rn^) < 


Cr, 


dmiwj’ 




( 8 ) 


I± = -j[p±{P± - a±) - p±{q± - a±)Y 
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where Oo± = H±(0). The invariants depend on the auxil¬ 
iary functions p± (scaling factors of the expansion modes) 
and a± (mass scaled centers of the dynamical modes of 
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the invariant). They satisfy the auxiliary (Ermakov and 
Newton) equations 


w + = (11) 

P± 

d± + = po±- (12) 


Dynamical expansion modes (not to be confused 

with normal modes) may be found. These are exact time- 
dependent solutions of the Schrodinger equation and also 
instantaneous eigenstates of the invariant Q, 
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where cr± = and $ra(cr±) are the eigenfunc¬ 

tions of the static harmonic oscillator at time t = 0 . 
Within the harmonic approximation the NM wave func¬ 
tions ji/ilO evolve independently with H". They may 
be written as combinations of the expansion modes, 
|' 0 (t:(O) = Sn Cn±|V'”±) with normalized constant ampli¬ 
tudes. The average energies of the n-th expansion mode 
for two NM are 
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In numerical examples the initial ground state is, in the 
harmonic approximation, of the form |'i/'o+(O))|' 0 o-(O)); 
so the time dependent energy is given by E"{t) = -|- 

Eq_ . Note that if we impose both unitary operators U(t) 
and U{t) to be 1 at t = 0 and tf, the transformed wave 
function \tjj") and the laboratory wave function lip) will 
be the same at both these times and the energy E"{t = 
0,tf) will be the same as the laboratory-frame energy. 
Both unitary transformations satisfy this provided that 
wi( 4 ) = 0 , where 4 = 0 ,t/, as long as the quadratic 
approximation in the Hamiltonian ([S]) is valid. 

For a single harmonic oscillator without the indepen¬ 
dent term in Eq. m, i.e. with a fixed center, the 
frequency in a trap expansion was already inverse en¬ 
gineered in [ 2 ^. For this case we can use the same no¬ 
tation as before but no subindices for the auxiliary func¬ 
tions. a is zero for all times, and in the Ermakov equa¬ 
tion the conditions p( 0 ) = 1 , p{tf) = 7 = and 

p{th) = pith) = 0 , suffice to avoid any excitation (since 
[H{tb), I{tb)] = 0 ) and ensure continuity of the oscilla¬ 
tor frequency. Any interpolated function p(t) satisfying 
these conditions provides a valid f2(t). Similarly, in har¬ 
monic transport of an ion (with the trap moving rigidly 
from 0 to c? with a constant frequency [7|) the auxiliary 
equation for p becomes trivially satisfied by p = 1 and, 
to avoid excitations and ensure continuity, a may be any 
interpolated function satisfying Q(( 0 ) = 0 , a(tf) = d, 



FIG. 1: (Color online) Final excitation energy at t/ for the 
expansion of two ^'’Ca"'" ions, with respect to the final ground 
state (quantum) or the final equilibrium energy (classical). 
The initial state is the ground state (quantum) or the equi¬ 
librium state (classical) for the initial trap. The dashed red 
line is the excitation in the harmonic approximation, using 
Eq. (1141) for the NM energies, with the protocol obtained 
by the shooting method; the solid blue line (classical) and 
black triangles (quantum) are for the same protocol but with 
the dynamics driven by the full Hamiltonian O- The dotted 
green line is for the protocol Gnu with the full Hamiltonian. 
The parameters used are loq /{ 2 ' k ) = 1.2 MHz and 7 ^ = 3. 


a{tb) = a{tb) = 00- Instead of these simpler settings, 
when inverse engineering the expansion of the ion chain 
the auxiliary equations m are non-trivially coupled and 
have to be solved consistently with Eq. m, since 
and po± are functions of the same frequency wi. In other 
words, only interpolated auxiliary functions p±(t), a±(t) 
consistent with the same uJi(t) are valid. 

For both NM, we impose for Eq. (ED the bound¬ 
ary conditions (BC) p±(0) = I, p±{tf) = 7 , p±{tb) = 
p±{tb) = 0. Here ujq = wi(0) and ujf = uJi{tf). The BC 
for the second set of equations are a±(tb) = a±{tb) = 
a±{tb) = 0. Eq. (IT^ with Eq. ([5|) implies that at the 

boundaries we must have | ~ ^i(if)) = 0. This is 

satisfied by imposing ioiitb) = 0, wi(tb) = 0. Substitut¬ 
ing these conditions in Eq. ED we finally get the extra 
BC pf(tb) = pfitb)=Q. 

To engineer the auxiliary functions we proceed as fol¬ 
lows: first we design p-(t)^ so as to satisfy the 10 BC 
for p-{tb) and their derivatives. They could be satisfied 
with a ninth-order polynomial, but we shall use higher or¬ 
der polynomials so that free parameters are left. These 
may be chosen to satisfy the equations for the remaining 
BC for a± and p+. uji{t) is deduced from the polyno¬ 
mial using Eq. ED SO it becomes a function of the free 


^ We choose p_ instead of p+ since $1+ > n_. The effective trap 
for the plus (+) mode is thus tighter and less prone to excitation 
than the minus (—) mode. Designing first p_ guarantees that 
this ‘weakest’ minus mode will not be excited. 
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FIG. 2: (Color online) The values of the optimizing free pa¬ 
rameters in the two-ion expansion aio (solid blue line) and an 
(dashed red line) in the expansion of two ^^Ca'*' ions starting 
in the ground state. ujq /{ 2 ' k ) = 1.2 MHz, 7 ^ = 3. 



tf (jws) 

FIG. 3: (Color online) Final excitation energy for the expan¬ 
sion of a ion chain (solid blue line) and a ^Be"*"- 

“^^Ca^-^Be^ chain (dashed red line) starting in the equilib¬ 
rium configuration. The protocols are optimized with four 
(for ®Be+-‘‘°Ca"'') and two (for ®Be^-'*°Ca'''-®Be''') free pa¬ 
rameters, see the main text. cJo/(27r) = 1.2 MHz, 7 ^ = 3. 


parameters. There are different ways to fix the free pa¬ 
rameters so as to satisfy the remaining BC and design 
the other auxiliary functions. In practice we have used 
a shooting method [s^. The BC used for the shooting 
are a±(0) = q;±( 0) = p+ = 0 and p-i-(O) = 1. Note 
that if Qf±(tb) = 0 , then a.±{tb) = 0 since we impose 
wi(ib) = wi(ib) = 0. The differential equations (ITT1) for 
p+{t) and ( 1121 ) for a± are now solved forward in time. 

In the following one must distinguish between single¬ 
species and mixed-species ion chains. A consequence of 
having equal mass ions is that a-{t) is 0 at all times (be¬ 
cause the ion chain is symmetric, and thus the center of 
mass remains static) so we only have to design the three 
auxiliary functions p±{t) and a+{t). When both ions are 
of different species, the chain is not symmetric anymore, 
so we also need to design a- taking into account its BC. 

The MatLab function ‘fminsearch’ [s^ is used to find 
the free parameters that minimize the total final energy 
for the approximate Hamiltonian, f) + f), see 

Eq. (fT4l) . For equal mass ions, an 11-th order poly¬ 
nomial p-{t) = X]n=o parame¬ 

ters, is enough to achieve negligible excitation in a range 
of times for which the harmonic approximation is valid. 
Only two free parameters are needed to satisfy the BC 
a+{tf) = a.+ {tf) = 0 , whereas p+(tf) = 7 is also nearly 
satisfied for all values of these free parameters because 
the evolution of this scaling factor is close to being adi¬ 
abatic. LUilt) is then a function of the free parameters 
oio, Oil. Fig. [IJdepicts the final excitation energy for op¬ 
timized parameters in the harmonic approximation, us¬ 
ing Eq. (|n]), and with the full Hamiltonian ([T]), whereas 
in Fig. [2] the values of the optimizing free parameters 
are represented. The quantum simulations (triangles in 
Fig. m are performed starting from the ground state of 
the Hamiltonian © at t = 0 , which is calculated nu¬ 
merically. For the corresponding classical simulations we 
solve Hamilton’s equations for the two ions in the lab¬ 
oratory frame with Eq. ©: the excitation energy is 


calculated as the total energy minus the minimal energy 
of the ions in equilibrium. The initial conditions corre¬ 
spond as well to the ions in equilibrium. As the potential 
is effectively nearly harmonic and the evolution of wave 
packet’s width {p±) is close to being adiabatic, the classi¬ 
cal excitation energy reproduces accurately the quantum 
excitation energy, as demonstrated in Fig. [TJ Quantum 
calculations are very demanding, in particular with three 
or more ions, so that we shall only perform classical cal¬ 
culations from now on. 


For two different ions, we use a 13-th order polyno¬ 
mial p-{t) = which is enough to nearly 

satisfy a±{tf) = a±{tf) = 0 and p+{tf) = 7 by finding 
suitable values for the four free parameters aio- 13 . As 
before, p+(tf) = 7 is nearly satisfied without any special 
design. Fig. [3] shows the final excitation for a chain of 
two different ions. The excitation is higher than for equal 
masses. Both for the equal mass and different mass ex¬ 
pansions, the (exact) excitation energy increases at short 
times, where the quadratic approximation to set the NM 
Hamiltonians fails, see Figs. [T]and[31 Further simulations 
indicate that the larger the ratio between the masses, the 
higher the excitation. 


A less accurate, approximate treatment is based on the 
simpler polynomial ansatz p_ = X]n=o without free 
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FIG. 4: (Color online) Final excitation vs. final times for ex¬ 
pansions of two equal ions (solid blue line), 4 equal ions (red 
dots) and 8 equal ions (dash-dotted black line). The simu¬ 
lations are performed according to the approximate protocol 
in Eq. d) and by solving the classical equations of mo¬ 
tion for '‘“Ca"*' ions. The initial ion chain is at equilibrium. 
ujo/{2'k) = 1.2 MHz, 7 ^ = 3. 


parameters,^ 

p- = 126(7-l)s^-420(7-l)s®-h 540(7-l)s^ 

- 315(7-l)s®-h 70(7-l)s®-hi, (15) 


s = t/tf. While the BC of p+ and a± are in general not 
accounted for exactly, an advantage of this procedure is 
that there is no need to perforin any numerical minimiza¬ 
tion. This is useful to generalize the method for larger 
ion chains. For equal masses, both a- = 0 and p-(t) 
are correctly designed, so that the center of mass is not 
excited. From Eq. m, oJi (t) is given by 



(16) 


where A- = is a constant, see Eq. ([3|). In Fig. [T] 

we compare the performance of this approximate proto¬ 
col and the one that satisfies all the BC in the two-equal- 
ion expansion. 


III. N-ION CHAIN EXPANSION 

We now proceed to extend the results in the previous 
section to larger ion chains governed by the Hamiltonian 


^ As in the transport of two ions [^, an alternative ansatz to 
the polynomial is p-(t) = 1^ + cos , 

where an = (-150,25,-3). In numerical calculations the poly¬ 
nomial ansatz CSli performs slightly better than the cosine-based 
one. 


FIG. 5: (Golor online) Classical trajectories of eight expand¬ 
ing '^°Ca+ ions. The evolution is performed according to Eq. 
(TTOl) . 6L)o/(27r) = 1.2 MHz, 7 ^ = 3, t/ = 4.4 ps. 


L The equilibrium positions can be written in the form 




where 


l\t) = 


Cc 

uo{t) 


and the Ui are the solutions of the system 

i-l , N 


u. 


-E f ^2 + E ^ 

^ {Ui - Ujf ^ 


{Ui - UjY 


= 0 . 


The NM coordinates are thus defined as [if 

Jo)^ 


= E a«v^(9* - ). 


(17) 


(18) 


(19) 


( 20 ) 


where the NM subscript v runs now from 1 to TV. 
Conventionally the v are ordered from the lowest to 
the highest frequency [s^. As for two ions we de¬ 
fine E(gi , g2, d3, <In) as the coordinate-dependent part 
of the Hamiltonian O- The a^i are the compo¬ 
nents of the ly-th eigenvector of the symmetric matrix 

that, together with the 

eigenvalues Xi, = will usually be determined nu¬ 
merically [ 3 ^. They are normalized as — 1- 

As uq is common to all ions, it can be shown that 
Hi,(t) = where A^ is a constant. 

Generalizing the steps leading to Eq. the Hamil¬ 
tonian in a NM frame up to quadratic terms becomes 


^' = E 


Pi/ I ^02^2 


POi/Pi. 


( 21 ) 


where the Pi/ are momenta conjugate to the q^, and 
Poi^ = - ■ As for two ions, all the poi. 
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FIG. 6 : (Color online) Comparison of final excitation qnanta (classical simulation as in Fig. 1) vs final time in the expansion 
of 2 ^°Ca+ ions following the shooting protocol (solid blue), linear protocol in Eq. (12411 (dotted green) and cosine protocol in 
Eq. (I25II (dashed red). In a) we plot in logarithmic scale up to times where only the shooting protocol reaches the level of 0.1 
excitation qnanta. In b) we extend the analysis up to longer final times, so that the best of the cosine protocol reaches also 0.1 
final excitation quanta. a;o/(27r) = 1.2 MHz, 7 ^ = 3. 


are proportional to . We now apply the unitary 

transformation U = and find the effective 

Hamiltonian 


H" 


E 





2 




( 22 ) 


This Hamiltonian is similar to the one for two ions 
The corresponding set of auxiliary equations is also sim¬ 
ilar to Eqs. (El) and m, 


Pv + ^tpv — 


+ — Poi/. 


(23) 


The BC for inverse engineering read Pv{0) = 1, Puitf) = 
7 , pv{tb) = Pv{tb) = 0 , av{h) = au{tb) = a„{tb) = 0 . 
When introducing the BC for the a^, in the set of New¬ 
ton’s equations, we get from all of them the same con¬ 
dition + (jJi{tb) = 0, which is satisfied for uJi{tb) = 

LOl{tb) = 0 . 

Fig. 0] depicts the excitation for expansions of single¬ 
species ion-chains, with approximate (non-optimized) 
protocols that use Eqs. ED and ED, but with the low¬ 
est frequency mode, ^ = 1, instead of the minus {—) 
mode. The longer the chain the lower the fidelity of the 
protocol, as more terms are neglected in the NM approx¬ 
imation and more boundary conditions are disregarded. 
However, the protocol still provides little excitation at 
long enough final times in the most demanding simula¬ 
tion that we examined, N = 8. Fig. [5]shows the position 
of the ions, and the trap frequency along the evolution 
time for the eight-ion chain, ending up with a separation 
between ions twice as large as the initial one, in times 
shorter than 4 ps (Fig. 0D without any significant final 
excitation. 

In Fig. [ 3 ] the excitation for an expansion of the 
two-species chain ®Be“''-"^°Ca+-®Be“'' is depicted. The 


minimization technique was used with two free param¬ 
eters, that is, with an 11 -th order polynomial ansatz for 
Pi,=i(t). The excitation is smaller than for the shorter 
chain ^Be^-^^^Ca^ (with a 13-th order polynomial for p_) 
due to the symmetry in the three-ion chain, which leaves 
two of the NM static and unexcited. 


IV. DISCUSSION 


We have designed fast diabatic protocols for the time 
dependence of the trap frequency that suppress the final 
excitation of different ion-chain expansions or compres¬ 
sions. Unlike the simpler single-ion expansion [ 2 ^, the 
inverse design problem of the trap frequency for an ion 
chain involves coupled Newton and Ermakov equations 
for each dynamical normal mode. We found ways to 
deal with this inverse problem by applying a shooting 
technique in the most accurate protocols, and effective, 
simplifying approximations. 

These protocols work for process times for which the 
quadratic approximation for the Hamiltonian is valid. 
Longer and more asymmetric chains need larger times 
than shorter and symmetrical ones. The examples show 
that these times are compatible with current quantum 
information protocols, so many processes may benefit by 
the described trap frequency time dependencies. 

The designed protocols provide a considerable im¬ 
provement in final time and excitation energy with re¬ 
spect to simple, naive protocols. For the expansion of 
two considered in Fig. 1 we compare in Fig. [ 6 ] 

the excitation energy of the shooting protocol with two 
simple protocols that drive the frequency uii linearly. 


U!i(t) — Wo + 


UJf-UJo, 

tf ’ 


(24) 
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and following a cosine function, 

CjQ+UJf UJQ-UJf /ttA 

The simulations are classical, as described in Sec. El 
Fig. [6] (a) compares the excitations at short times. For 
tf ^ 2.5 fis, the shooting protocol reaches a low excita¬ 
tion of 0.1 vibrational quanta, four orders of magnitude 
smaller than the excitations due to the simple methods. 
In Fig. [6] (b) the excitations are represented for longer 
protocol times. The smoother cosine protocol behaves 
better than the linear one and finally reaches an excita¬ 
tion of approximately 0.1 quanta for ty ~ 20 fis. 
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